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We discuss reorganizing finite temperature perturbation theory using 
hard-thermal-loop (HTL) perturbation theory in order to improve the con- 
vergence of successive perturbative approximations to the free energy of a 
gauge theory. We briefly review the history of the technique and present 
new results for the three-loop HTL-improved approximation for the free 
energy of QED. We show that the hard-thermal-loop perturbation reor- 
ganization improves the convergence of the successive approximations to 
the QED free energy at intermediate coupling, e ~ 2. The reorganization 
is gauge invariant by construction, and due to cancellation among various 
contributions, one can obtain a completely analytic result for the resummed 
thermodynamic potential at three loops. 

PACS numbers: ll.15.Bt, 04.25.Nx, ll.10.Wx, 12.38.Mh 

1. Introduction 

In the early 1990s the free energy of a massless scalar field theory was 
calculated to order g 4 in Refs. [lj [2]. This was quickly followed by similar 
calculations in QED [3j and QCD [2]. The scalar, QED, and QCD free 
energies to order g 5 were then obtained in Refs. [U [5], Refs. [6j [7] and 
Refs. [3 [9], respectively. Recent results have extended the calculation of the 
QCD free energy by determining the coefficient of the g & log(<?) contribution 
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Fig. 1. Successive perturbative approximations to the QED pressure (negative of 
the free energy). Each band corresponds to a truncated weak-coupling expansion 
accurate to order e 2 , e 3 , e 4 , and e 5 , respectively. Shaded bands correspond to 
variation of the renormalization scale \x between ttT and A-kT . 

|10j . For massless scalar </> 4 the perturbative free energy is now known to 
order g 6 p] and g 8 log (5) [12]. 

However, the resulting weak-coupling approximations, truncated order- 
by-order in the coupling constant, are poorly convergent unless the coupling 
constant is extremely small. For example, simply comparing the magnitude 
of low-order contributions to the Nf = 3 QCD free energy one finds that 
the <7g contribution is smaller than the g 2 s contribution only for g s < 0.9 
(a s < 0.07). This is a troubling situation since at phenomenologically acces- 
sible temperatures near the critical temperature for the QCD deconfinement 
phase transition, the strong coupling constant is on the order of g s ~ 2. 

The poor convergence of finite-temperature perturbative expansions of 
the free energy is not limited to QCD. The same behavior can be seen in 
weak-coupling expansions in scalar field theory \Y'i\ [T5] and QED [6]. In 
Fig. Q] we show the successive perturbative approximations to the QED free 
energy. As can be seen from this figure, at couplings larger than e ~ 1 the 
QED weak-coupling approximations also exhibit poor convergence. For this 
reason a concerted effort has been put forth to find a reorganization of finite- 
temperature perturbation theory which converges at phenomenologically 
relevant couplings. 
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There are several ways of systematically reorganizing the perturbative 
expansion to improve its convergence and the various approaches have been 
reviewed in Refs. |15[ [TBI E]- Here we will describe recent advances in 
the application of hard-thermal-loop perturbation theory (HTLpt) [HI [TjJl 
[201 l2Tj I22j . The HTLpt method is inspired by variational perturbation 
theory [231 [Ml EHJ [26] and is a a gauge- invariant extension of scalar screened 
perturbation theory (SPT) [22 [281 El E] . The basic idea of the technique is 
to add and subtract an effective mass term from the bare Lagrangian and to 
associate the added piece with the free Lagrangian and the subtracted piece 
with the interactions. However, in gauge theories, one cannot simply add 
and subtract a local mass term since this would violate gauge invariance. 
Instead one adds and subtracts a HTL improvement term which modifies 
the propagators and vertices in such a way that the framework is manifestly 
gauge-invariant. The free part of the Lagrangian then includes the HTL 
self-energies and the remaining terms are treated as perturbations. 

In this brief proceedings review we present results of a calculations of the 
QED free energy (pressure) to three-loop order in HTLpt based on the work 
detailed in Ref. [29]. As we will show, the next-to-leading order (NLO) and 
next-to-next-to-leading order (NNLO) HTLpt resummed QED free energy 
give approximations which show improved convergence for couplings as large 
as e ~ 2.5 (see Fig. [2]). In addition, we compare our results to those obtained 
using the 2PI ^-derivable approach [30] and show that at three loops the 
agreement between the HTLpt and ^-derivable approaches is quite good. 



2. Formalism 

The Lagrangian density for massless QED in Minkowski space is 

£qed = --F^F^ + iftfDrf 

+£ g { + £ g h + A£qed . (2.1) 

Here the field strength is F^ v = d ii A v — d v A^ and the covariant derivative 
is = + ieA^ . The ghost term £ g h depends on the gauge-fixing term 
£ g f. We will use dimensional regularization with a renormalization scale \i 
and covariant gauge fixing such that the ghost terms decouple. 

Hard-thermal-loop perturbation theory is a reorganization of the pertur- 
bation series for thermal gauge theories. In the case of QED, the Lagrangian 
density is written as 



£ = (£qed + £htl) 



+ A£ H tl 

Voe 



(2.2) 
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The HTL improvement term is 




+(l-5)im}^{-^) c. (2.3) 

where = (1, y) is a light-like four-vector, and (. . .) y represents an aver- 
age over the directions of y. The term (|2.3h has the form of the effective 
Lagrangian that would be induced by a rotationally-invariant ensemble of 
charged sources with infinitely high momentum. The parameter mp can 
be identified with the Debye screening mass and the parameter mj can be 
identified as the induced finite-temperature electron mass. HTLpt is defined 
by treating 5 as a formal expansion parameter and expanding order by or- 
der in 5 around 5 = 1. This generates loops with fully dressed propagators 
and vertices and also automatically generates the counterterms necessary 
to remove the dressing as one proceeds to higher loop orders. 

If the expansion in S could be calculated to all orders, the final result 
would not depend on mr> or mj. However, any truncation of the expansion 
in 5 produces results that depend on mp and mj. Some prescription is 
required to determine and mj as a function of T and e. For example, 
one can choose to treat both as variational parameters that should be de- 
termined by minimizing the free energy or one can fix mc and mf using a 
perturbative prescription. We will compare both methods. We will obtain 
the thermodynamic potential £l(T, e, mo, mj, fj,, 5 = 1) which is a function 
of the mass parameters mo and mf. The free energy T is obtained by 
evaluating the thermodynamic potential at the appropriate values of the 
thermal masses. Other thermodynamic functions can then be obtained by 
taking appropriate derivatives of T with respect to T. 

3. Results 

In this section we present the final renormalized thermodynamic poten- 
tial explicitly through order 5 2 , aka NNLO, as obtained by us in Ref |29j . 
The final NNLO expression is completely analytic; however, there are some 
numerically determined constants which remain in the final expressions at 
NLO. 

3.1. Next-to-leading order 
The renormalized NLO thermodynamic potential is 

nlo = ll + -N f -15m 3 D -- h og | --+ 7 + — I ml, 
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+60N f (vr 2 - 6) mj + N f - 
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+ 15mn 



+5 ^log | - 2.33452 j m 2 D - 45 Hog | + 2.19581 ) mj 
-30 Hog | - - + 7 + 2 log 2 W, + 180m D m 2 




(3.1) 



where we have introduced the dimensionless parameters rri£> = tiid/ (2ttT), 
rhf = rrif/(2irT) , and (i = ^/(2-kT). 



3.2. Next-to-next-to-leading order 
The resulting NNLO thermodynamic potential is 

7T 2 T 4 
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5.5. Free Energy 

The mass parameters mo and mj in hard-thermal- loop perturbation 
theory are in principle completely arbitrary. To complete a calculation, it 
is necessary to specify m_n and mj as functions of e and T. In Ref. [29] 
we considered two possible mass prescriptions in order to see how much 
the results vary given the two different assumptions. First we considered 
the variational solutions for the thermal masses and second we considered 
using the e 5 perturbative expansion of the Debye mass [3TJ [7] and the e 3 
perturbative expansion of the fermion mass [32]. The resulting predictions 
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Fig. 2. A comparison of the renormalization scale variations between NLO and 
NNLO HTLpt predictions for the free energy of QED with Nf = 1 and the vari- 
ational Debye mass (left) and using the perturbative thermal masses (right). The 
bands correspond to varying the renormalization scale \i by a factor of 2 around 
fi = 2nT. 

for the free energy are shown in Fig. [2j As can be seen from these fig- 
ures both the variational and perturbative mass prescriptions seem to be 
consistent when going from NLO to NNLO. As a further check of our re- 
sults in Fig. we show a comparison of our NNLO HTLpt results with a 
three-loop calculation obtained previously using a truncated three-loop $- 
derivable approximation [30] . As can be seen from this figure, there is very 
good agreement between the NNLO <I>-derivable and HTLpt approaches out 
to large coupling. 

3.4- Conclusions and Outlook 

In this paper we discussed reorganizing finite temperature perturbation 
theory using HTLpt in order to improve the convergence of successive per- 
turbative approximations to the free energy of QED. We presented results 
of a recent three-loop HTLpt calculation of the pressure in QED and 
showed that the HTLpt reorganization improves the convergence of the suc- 
cessive approximations to the QED free energy at intermediate coupling, 
e ~ 2. We studied two different mass prescriptions and showed that the 
results for the free energy using both prescriptions were the same to an ac- 
curacy of 0.6% at e = 2.4. We also compared the HTLpt three-loop result 
with a three-loop ^-derivable approach [30] and found agreement at the 
subpercentage level. In closing, we mention that the HTLpt reorganization 
is gauge invariant by construction we were able to obtain a completely ana- 
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Fig. 3. A comparison of the predictions for the free energy of QED with Nf = 1 
between three-loop "^-derivable approximation [30] and NNLO HTLpt at /z = 2irT. 

lytic result for the resummed QED thermodynamic potential at three loops. 
This gives us confidence to apply the method also to full QCD. 
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